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ABSTRACT 


This thesis presents two methods of solving the n-job m-machine 
job shop scheduling problem. The criterion for optimality is the 
minimization of the total time to process all jobs on all machines. 
The technological ordering of machines for each job is fixed, known, 
and nonrandom. 

The first method presented, a graphical method, indicates a lower 
bound and an upper bound on the optimal time to process all jobs on 
all machines. 

The second method is a branch and bound algorithm. In principle 
an optimal solution can always be determined by this method. Only 
limited computational experience is presented for the algorithm but 


some methods for efficient computation are suggested. 
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I. INTRODUCTION 


The problem considered is that of processing n jobs on m machines 
so that the time to process all jobs on all machines is a aii: 
Each job is processed in turn on each machine according to a fixed 
machine ordering. The machine orderings for the jobs need not be the 
same. 

The assumptions upon which solutions to the problem are based are 
taken from Refs. 1 and 2. They are: 

1. The time to process each job on each machine is fixed, known, 
and nonrandom. 

2. A job may not be processed by more than one machine at a time 
and a machine may not process more than one job at a time. 

3. Once processing of a job by a machine has begun, it may not be 
interrupted cae processing of that job is complete. 

Two approaches to the problem were considered. The problem was 
first represented geometrically. An attempt was then made to solve the 
problem graphically. The graphical method did not always provide an 
optimal solution but did provide an upper and lower bound on the time 
to complete all jobs. 

The second approach to the problem was a branch and bound method 
of solution. In principle the branch and bound method always yields 
the optimal solution to the problem. Computational efficiency of the 
branch and bound algorithm presented has not been examined in detail, 


but some methods for efficient computation are suggested. 





IIT. NATURE OF THE PROBLEM 


A. DEFINITION OF A SCHEDULE 

The properties of a schedule must first be defined. A nachane 
Schedule can be represented by a machine number (or letter) followed 
by an ordered n-tuple specifying the job ordering on that machine. 
Since the processing time of each job on each machine is known, m of 
these representations together with the known machine order for each 
job constitute a complete schedule. For example, if jobs 1, 2, and 3 
are processed in the order 2, 3, 1 on machine A, then the representation 


ome ( 23,1). 


B. GRAPHICAL PRESENTATION OF A SCHEDULE 

A complete schedule can be represented on a graph called a Gantt 
chart [Ref. 3]: A Gantt chart is simply a means of presenting a 
schedule visually. A list of machines is presented as the vertical 
edge of the chart. The horizontal edge represents time. Beside each 
machine is a plot of the start and finish of each job on the time axis. 


An example of a Gantt chart is shown in figure 1. 


Machine 





Time 





In the example the representations A:(1,2,3), B:(3,1,2), and 
C:(1,3,2) are presented graphically. The machine orderings for the 
jobs and the time to process each job on each machine are also presented 
in the Gantt chart. In figure 1, job 1 must be processed esc aOm 
machine A for 6 units of time, then on machine C for 2 units, and 
finally on machine B for 5 units of time. The machine orderings and 
process times for jobs 2 and 3 may be obtained in the same manner from 


migure i. 


C. CRITERION FOR AN OPTIMAL SCHEDULE 
The criterion selected for determining an optimal schedule was to 
minimize the total time to process all jobs on all machines. 


The time to process any job 71 on all machines can be represented by 


m 
oy (cere) ct; Ss. 3 re 
j=1 ij 1 
where chee processing time for job 7 on machine j 
34S total idle time of job ij 
m = the number of machines 
and n = the number of jobs. 


The problem is then to minimize the time to process all jobs on all 


machines or, equivalently: 


Minimize [ Max { 


: a 2 ee as oe eae 


] VJ 1 


tums 


J 


It is necessary to preserve the ordering of machines for each job, to 
ensure that no two jobs are processed at the same time on any machine, 


and to prohibit any two machines from simultaneously processing any job. 





YD. PREVIOUS SOLUTION METHODS 

Many attempts have been made to develop a closed-form solution to 
the general n-job m-machine flowshop problem. Most of the solution 
methods devised are acceptable for small m and n but become computa- 
tionally inefficient when applied to large scheduling problems. 

A graphical method of solution of the 2-job m-machine problem has 
been proposed by Akers and Friedman [Ref. 4] and Szwarc [Ref. 5]. The 
graphical algorithm presented in these two papers was proved by Hardgrave 
and Nemhauser [Ref. 1] to provide the sequence of jobs for each machine 
which minimized the total time to complete both jobs. 

Brooks and White [Ref. 6] and Ignall and Schrage [Ref. 7] have 
investigated the branch and bound method in solving the flowshop prob!em. 
Ignall and Schrage's algorithm was applied to the case where the criterion 
for optimality was the minimization of the time to complete all jobs but 
the machine ordering for each job was identical. Brooks and White have 
Stated that their procedure uses a great deal of computation time ona 
computer. No detailed computational experience was presented for their 
problem. 

Integer programming has also been attempted as a method of solving 
the jobshop problem but, according to Conway, Maxwell, and Miller [Ref.3], 
very little progress has been made. Tne reasons for lack of success in 
this area appear to be the large size of the integer programs involved 
and the large time required to solve the resultant integer programming 


problem. 





ITI. GRAPHICAL SOLUTION TECHNIQUES 


A. HARDGRAVE-NEMHAUSER GRAPHICAL ALGORITHM 
1. Two-Job m-Machine Scheduling Problem 

Although Akers and Friedman [Ref. 4] and Szwarc [Ref. 5] have 
described the graphical method of solution, Hardgrave and Nemhauser 
[Ref. 1] have provided a rigorous proof that the algorithm leads to an 
optimal solution. For this reason, the graphical method is referred to 
as the Hardgrave-Nemhauser graphical algorithm. 

a. The Geometric Model 

A geometric model of the 2-job m-machine problem is used in 
the Hardgrave-Nemhauser graphical algorithm. A coordinate axis is 
drawn with the ordinate representing the time to complete one job on 
all machines, the abscissa the time to complete the other. The known 
machine orderings and their associated processing times are marked along 
the axis corresponding to the appropriate job. Lines are drawn perpen- 
dicular to the axes at the finish times of the machine processing 
operations. The rectangles formed by the intersection of the spaces 
corresponding to a particular machine for each of the two jobs are shaded 
imetas in Tiaureme) . 

From any point a line can be drawn parallel to or at a 45 
degree angle with the axes. A line of this type from (0,0) to (S)5S,) 
that does not pass through any infeasible region (shaded area) is called 
a feasible path. A feasible path specifies an ordering of jobs on 
machines. That is, if the path passes beneath an infeasible region, then 


the job on the ordinate is processed before the job on the abscissa on 





ORDER OF OPERATIONS (PROCESSING TIME) 
SOB Atz), B(4), €(3), DU} 
NOpec> D3), AlZ), C13). B(2) 
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the machine generating that infeasible region. Then, a feasible path 
from (0,0) to (Sy 5So) specifies a complete schedule. 

Hardgrave and Nemhauser proved that the shortest feasible 
path from (0,0) to (S, So) yields the schedule which minimizes total 
process time of all jobs on all machines. Note that travel for one unit 
of time in any direction indicates the manner in which jobs are processed. 
Travel in a vertical direction indicates that only job 2 is being 
processed, ina horizontal direction that only job 1 is being processed, 
and diagonal travel indicates simultaneous processing of both jobs. The 
length of the shortest feasible path, hy from (0,0) to (S1 4S) was 
proved to be: 


A = Sim * Pam * MaxtSy - Sin - Pams So - Sam) 
where es set up time on machine m for job 1, 1=1,2 
oa processing time on machine m for job i, i=1,2 
and m is the machine generating the infeasible region that 


causes a conflict between jobs 1 and 2. 
The equation for shortest path is based on the assumption 
that job 1 precedes job 2 on machine m. If job 2 precedes job 1, then 


the equation for the shortest path is 


Ay = Som + Pom * Max (S, — 50m 7 Pom? 39 ? Som) 


and the optimal ordering of jobs on machine m can be determined from 


oe Min(A) 5A). 


That is, if A, <2. then job 1 precedes job 2 on machine m. If 


] 2 


12 





aes dy then job 2 precedes job l onm. Finally if A then 


2 lee 
either can occur and there is an alternate optimal solution to the 
problem. 

b. The Graphical Algorithm 

The graphical algorithm of Hardgrave and Nemhauser consists 

of rules for branching the path around infeasible regions. The algorithm 
USP 
Stier |: Starting at (0,0) move diagonally until an infeasible region 


is encountered, say on its bottom edge at point (t, s (If the left 


om) 
side of the infeasible region is encountered, the rules are similar.) 
Sir 2: Branch in two directions around the infeasible region: 
(a) Horizontally along the bottom edge until the point (Ss), ae ae Son) 
is reached. The return to STEP 1, using (Sy, + Pips So) as the new 
Starting point. (b) From point (t, Som) move in reverse along the path 
). From this point move vertically along the 


to point (s, a 


m?>2m Tm 


left edge of the infeasible region until point (s)_> So + Pom) is 
reached. The return to STEP 1 using this as a starting point. 
SEER SO! When the top or right edge of the outer rectangle is en- 
countered, move along that edge to the finish point (S),S,). The 
shortest path is the smallest of those paths generated in this manner. 
From this path the optimal schedule is found. 
2. The n-Job m-Machine Problem 

The geometric model and the graphical algorithm were extended to 
the general case by Hardgrave and Nemhauser, The model became a graph 
with n axes, one for each job. The infeasible regions in the model were 


n-dimensional polyhedrons. An optimal solution was characterized by the 
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shortest line which passed from (0,...,0) to (S)5...55)) without passing 
through any of the infeasible regions. Branching rules were given in 
the paper. Even for n=3 , this procedure was difficult to actually 


perform. 


B. GRAPHICAL BOUNDS ON OPTIMAL PATH LENGTH 
1. Description of the Metnod 

The (5) 2-dimensional orthogonal projections of the n- 
dimensional polyhedrons are simply rectangles. The projections of the 
entire n-dimensional graph are \> | 2-dimensional graphs. If these 
2-dimensional graphs are looked upon as the graphs associated with 
Bb 2-job m-machine problems, the shortest path in each graph may be 
found in the manner described previously. Each of these shortest paths 
indicates the ordering of two jobs on each machine. Considering all of 
the 2-job solutions, there are a pairwise orderings of all the jobs 
on all the machines. These pairwise orderings of jobs on machines indi- 
cate an ordering of all jobs on all machines. This ordering of jobs on 
each machine may be the plotted as a schedule on a Gantt chart as in 
figure 1, since the other required information is known (i.e., job's 
machine ordering and processing times). The time to complete the sched- 
ule is then determined by completion of the Gantt chart. 

W. Szwarc [Ref. 2] has also arrived at this technique and called 
it an approximate solution to the n-job m-machine scheduling problem. 
He showed that, under certain conditions, this method yielded an optimal 
solution to the problem in which each of the jobs had the same machine 
orderings. 

The method described does not always guarantee a feasible ordering 


of jobs on machines. Likewise, if the ordering of all jobs on all 
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machines obtained is feasible, there is no guarantee that the ordering 
is optimal. In other words, the time to completion indicated by the 
Gantt chart is not necessarily the minimum time to complete all jobs on 
all machines. Of course the solution obtained in this way is sometimes 
optimal. 
2 ASSES en ay 

The ordering of jobs on a machine A is infeasible if, for any 
of the three jobs i, j, and k, the pairwise orderings take the form 
Pied), Atl jsk), Al(k,i). 

If the ordering of jobs on some machine is not feasible, the 
ordering may be made feasible by applying the following branching rule: 

For one machine suppose the ordering of jobs is infeasible. 
That is, for some machine A the pairwise orderings for three jobs i, j, 
and k are of the form A:(i,j), A:(j,k), A:(k,i). Then, in each of the 
three associated graphs, branch the path in each of the ah possible 
ways around the rectangle causing the infeasibility. That is, generate 


the orderings 


CRA ie epee ekls As Niek) 
(2) ASVisd)6 As(gk), A:(k,7) 
Sy REGie ghee Uksd)) 2A (15K) 
(4) A:(i,j), Az(k,j), A:(k,i) 
(Oya ee ede is As isk) 
CCN ieee leis Ackksa) 
Walmer (351). A(koa), As(isk) 
and (8) A:(j,i), Az(k,j), Az(k,i) 


Of these orderings, (2) is the original infeasible ordering and (7) is 


also an infeasible ordering. Therefore, these orderings are not 
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ponendencce First, choose those orderings that do not cause any of the 
other orderings to become infeasible. Next, among those orderings, 
choose the one which causes the least increase in each of the shortest 
paths in the three associated graphs. Although nothing can be said re- 
garding the optimality of this generated path, it is feasible. Certainly 
if n is large this method is not practical. 
3. Upper and Lower Bounds 

The time to complete the feasible schedule generated by the 
solution of the 5 | 2-job m-machine problems clearly provides an upper 
bound on the minimum time to process all jobs on all machines. Call 
this time M A lower bound can also be found, thereby providing an 
upper and a lower bound on the optimal solution. 
Let he = the shortest path for the 2-job m-machine problem with 


J 
coordinate axis labeled i and j. 


and A* = the minimum time to process all jobs on all machines. 
Then, * 2G ae for each 1,j), if) 
or tomers Max (ure). 49> = n 
— 1ifj 1j 5 5 RZ e@ ee 9 e 
Then, us = Max (A..) is a lower bound on A*, or 
Wd 2 ale 
dy SSeS di 


Clearly, if dy = dy then A* = i = dy and the optimal schedule has 

been generated. In most cases, however, this will not be the case and 

the method merely indicates an upper and a lower bound on the solution. 
Whether the upper bound is a tight bound or not, it can be used 


to eliminate searching some of the branches in the branch and bound 


method of section IV. 





4. An Example 


As an illustration of the method, consider the 3-job 4-machine 


problem with the machine orderings: 


Job Machine (process time) 
1 Pie emeco (7), DiG) 
2 BREW (07s) DAZ). Ce) 
and i) A(4), Ci2) 5 Die) B(4). 


The problem is to minimize the time to process jobs 1, 2, and 3 on 

: SF) eon... 
machines A, B, C, and D. Since n=3, there are (> = Saar ae 3 
2-job 4-machine problems to be solved graphically. The solutions to 
the three problems are shown in figure 3. The job orderings obtained 
from these three graphical problems are: 


Mie Jeb bys. Joumse\..=22, As(le3 


ls 
Be(133 
C35 | 
= 


(3 
(3 
iA) Salo 2 OS a8 viet dog=l5s es 
(2 
(3, 


DuUc.o 
CeelonG 
Bi(3 2 


) 
) 
2) 
) 
) 
) 


(3) Job 1 vs. Job 2, 721, Ae(1;2 
eee |} 
Celie) 
DRCAaIe ye 

The aggregate job orderings formed from these three pairwise orderings 


ones 





job 3 


job 3 





machine 





S) 


10 


figure 3 
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RG) 
Becta. ) 
C213 22) 
and eecracral lis 


These orderings are all feasible. The Gantt chart is drawn and the 
upper bound on the time to process all jobs is obtained (A =22). The 
lower bound is dy =Max( y 32 An32 Ayo) #22. Since Mah =22, the solution 
4*=22 1s optimal and the optimal schedule is presented in the Gantt 


Giart in figure 3. 
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IV. BRANCH AND BOUND SOLUTION 


A, PROBLEM FORMULATION 

Consider again the problem of processing n jobs on m machines so 
that the time to process all jobs on all machines is a minimum. The 
machine ordering for each job is fixed and is known, as is also the 
processing time for each job on each machine. 

As stated previously, the total time to process any job i on all 


machines can be represented by 


m 
Pa (t; 5) + Sey 7 = Ty..egn 
where tis = the processing time for job i on machine j, 
is total idle time of job i, 
m = the number of machines, 
and n = the number of jobs. 


The problem is then to minimize the time to process all jobs on all 
machines. This can be written as 


m 
Minimize [ Max ope yee 


: 3) at Si, 1 = 1,....n } ] 
j=] : 


Uy 

At the same time it is required that the ordering of machines for each 

job be preserved, no two jobs be processed at the same time on any 

Machine, and no job be simultaneously processed by any two machines. 
The total processing time for job i can also be represented as 


m-| 


Be pet) eS 
fe] iF j im 
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where m is the number of the last machine, in order, to process job i 


and the other variables are as defined previously. Then, setting 


m- | 
Se = : get) Mie ae Sasha 
m j=1 1) 1 


a new form of the objective function is 


Minimize [ te (OS ime: mine ise enn pd 


aa is called the earliest start time of job i on machine m, or, the 
earliest start time of the last process which job i must undergo. 

Let Fey be the latest completion time of job i on machine 1. Since 
job i and job j cannot be processed at the same time on any machine 1, 
then 


sy |e 


4 


ail 
if job j precedes job 1 on machine 1. 


Additionally, if job i is to be processed on machine 1 before beijing 


processed on machine k, then 


Sik 2 Fay 
The original problem can now be restated as 


Minimize [ MX {s., +t. , i= 1,..n}]. 


m 


SUbVeeE LO > 


|v 


“1 Fare if j precedes i on |] 


or 33) 


lv 


Faas iieimenecedes. jvom linn tor all 1. Jaloe..5n 


eens 


pore ln 


21 





and ee oe ae 
1515 at 151] 


se Swe Ibe oe 
Ul —— il 


—N 
|v 
=— 


- fOr valid 4 


where the machine ordering specified for job i is TyoeeeoT,. 


B. BRANCH AND BOUND ALGORITHM 
1. Branching Method 

The problem, as stated, can be solved using a branch and bound 
algorithm similar to that used by Little and others [Ref. 8] to solve 
the traveling salesman problem. In the job shop scheduling problem a 
branching operation can be defined through the ordering of jobs on 
machines. First consider job 1 to be scheduled on all machines. Then 
define a branch by ordering job 2 and job 1 on the first machine to 
process job 2. These orderings would be 2< ] and 1< 2. Then 
calculate a lower bound on the time to process all jobs on all machines 


(i.e., Max (Si + te 
m m 


smallest lower bound, branch again in the same manner for the next 


, i=l,...,n) ). From the node giving the 


machine to process job 2. Continue branching in this way until both 
jobs 1 and 2 are scheduled on all machines. When this has been done, 
branch in three ways from the node with the smallest lower bound. That 
iemondnch on tnpesonderings o <1 <2, 1<3 <2 ,and | <2<43 

for the first machine to process job 3. These are the orderings if job | 


preceded job 2 on this machine. If on the contrary, job 2 preceded job | 


Ze 





on the machine in question, then the orderings would be 3 <2 <1, 
2<3 <1 ,and 2 <1 <3. Then calculate the lower bound and continue 
to branch on the smallest existing lower bound. Continue to add jobs 
until all jobs have been processed on all machines. When the smallest of 
the lower bounds has been found and all jobs have been processed on all 
machines, the optimal solution has been found. 

An illustration of the branching method is shown in figure 4. 
In the illustration, the set of all job-machine orderings is partitioned 
by ordering job | and job 2 on machine A. This is represented by 2 < ] 
and 1 <2. The fixed ordering of machines specified for job number 2 
is used for the two job partitioning. In the illustration, suppose that 
the ordering of machines A, B, C is the ordering specified for job 2. 
Job 1 must also be processed according to its specified ordering on 
machines. This is done by considering job | to be scheduled already. 
The ordering of jobs may be the same or different. By branching on 
successive orderings of jobs on machines, this technique is capable of 
representing any finite number of jobs and machines. It is limited by 
the number of computations. 

2. Lower Bound Computation and Algorithm 
If a job j follows a job 1 on machine 1, then it must be true 


that 


Likewise, if job i is processed first on machine k and then on machine I, 


then 


ts 
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py IEE Machine 


Level 


2< 1) 1<2 A 


2 az B 


Cie Ona) Ges) » 


BRANCHING METHOD 
ORDERING OF JOBS ON MACHINES 


figure 4 
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It is possible to record these times in tableau form and, consequently, 
to know the time required to process all jobs on all machines (through 


the present branch). This time is S; Baw ct t. a where the 1s the last 
mn om : 


machine encountered by job i. 


The tableau to be used is one of the form illustrated in figure 


151 Vg 1 


oeevniere 1. - S. . +t. ., 122 1,...,n and S.. and F.. are as previ- 
' m m ws "J 


ously defined. 





The methods of computing the required lower bound and describing 
the branching operations is given in the following algorithm. 
Ser <1: For the initial tableau, fill in the blocks of the tableau 
as if each job were to be processed independently of every other job. 


When completed, every row i should have entries of the form 


2 ye eget i eos 


and F. . =T7 Tate, Wi ecceen ets 
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This schedule will undoubtedly be infeasible for joint processing 
of all the jobs for any interesting problem. If this schedule is feasible, 


then it is optimal and the minimum time to process all jobs is 


Max le. 
‘ { T., tem aleereees His Fy 


If this schedule is infeasible, go to STEP 2. 
See 2: For the first partition (suppose it is i < j , or i precedes 


j on machine k) evaluate the lower bound on 


Max 


It is known from the constraints that for job i to precede job j on 


machine k, it is necessary that 


The procedure, therefore, is: 


(a) If ene Says add ee Sip to S ik and to every S51 where ] 
follows k in the machine ordering for job j. That is, enter in the next 
tableau: 
34) - Fey - Sak a 351° for all 1 where 34] 2 Sax 
lee 5] Sees 
J Ue Jeaies Jed 

and Fay = 33) = tay> for all 1 where 33) 2 Sa: 


° { lites i lees, <A } ° 
1 1 


Gos 1? Fey SOE make no change. Record the lower bound, be {T.} : 
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Ser 3 When the lower bound for the branch i< j has been found, 
compute the lower bound for the branch j < i. in the same manner. 

STEP 4: Compare the lower bounds thus obtained and continue to branch 
along the path with the lowest computed lower bound. At each iteration, 
compare the latest finish time of the first job in the ordering with the 
earliest start time of the second job in the ordering. Then, if any 
change is necessary, compute the new Si Pi and and compare that 


Fay with the next job's Spy in the precedence ordering for that machine. 


Continue this procedure until, for the machine level reached, all 
ee. Faye lt Jcion., | 


Or S., > Fea if i<j on ] 


jl 
fOgmarieier) = l4.05,n 
and ] = Tyoeeesdy, 
and Se 
VO UG) = VsVy 
Sie ew _— 
De Um UALR Olea We listens Il, 


Record the lower bound, or {T.3 


SINE aye Continue to branch along the paths for which the least lower 
bound has been recorded until all jobs have been ordered on all machines. 
Then if there is no recorded lower bound in any branch that is less than 
that finally computed, the optimal solution is specified. 


3. An Example 


Consider again the example problem of scheduling three jobs, 1, 2; 


2] 





and 3, on four machines, A, B, C, and D. The machine orderings and 


processing times are: 


Job Machine (process time) 
1 A(3), B(3), C(7), D{6) 
2 B(5), A(6), D(2), C(2) 
2 PNA C2) eS) sf iB (4) 


The first tableau is 





On the first branch, since job 2's machine ordering is B, A, D, C, 
the set of all orderings is partitioned into 2< 1 and 1< 2 0n 


machine B. Illustrated, this is 
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Suppose the branch 2 <1 is checked first. As described in 


the algorithm, F op is compared with 1B" Note that F op =o | abe 1S. 


Fon > 518: Therefore a change in the tableau is made of the form: 


“Aa 


toe SS ee be 


J Jodp, dad 


Psy = Sie 


and jl 


are 1 Or ahilewnerce  S.- > Sa.% 
jl jl —~gk 


Or, since B is the second machine in job 1's ordering, 1 = B,C,D. 


Therefore ns 
“1B” fan “1p * 1B? 
aes: al alle 
Som = lo = iy = Sp was 
T= Sqp t hp = 2 
ig yay Fine < 
Ae = Sic t Hen 1 

and a a + tip = 2] 


The new tableau is then 


Max ¢ 7 3 
7 1 


(2) 2 
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At this iteration _ ie = 21. Record this as the lower bound 
and check the branch 1 <2. The same procedure is followed and the 


tableau is 


Max ¢ 7 3} 
; i 





Since both computed lower bounds are the same, the next branch 
can be made from either. Suppose 1 <2 on .B is chosen and partitioned 
again into 2 <1 and 1 <2 onA. Note that job 2 cannot precede 
job 1 on machine A with the specification that job 1 precede job 2 on 
machine B. That branch can, therefore, be omitted. A lower bound of 
o is then assigned. In the example, branch 1 <2 on A is then made. 
Noting that Fin < Sons no change is made in the tableau, but another 
branch is made, partitioning on 2 <1 and 1 <2 onD. This time, 


Owes ls > 4p and a tableau change is made. The new 


oD 
oe { LP } = 25, Likewise, for the branch 1 <2 onD, it is found 
that ia {T,} = 23. Since both of these bounds exceed 21 (the bound 
on 2 <1 on A), that bound determines the next branch. It was not 


necessary to compute other than Ts Sco Of TY, = 23. If all the bounds 


computed exceed 1). the new tableau could then be reconstructed. 
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Branching on 2< 1 on A, it is noted that T, = 32, so the 


branch 1 < 2 on A can be checked. It is noted that F S 


ein, 2° 
the same tableau is applicable for the next branch and the lower bound 
is still 21. The computations are continued in this manner until the 
two jobs are scheduled on all machines. Then machine 3 must be considered. 
The computational method is the same except that there are now three jobs 
to be considered. The branching and bounding is continued until the 
branch 2< 1< 3 onB is reached and its lower bound is computed. Since 


this lower bound is 22, a number smaller than all other computed lower 


bounds, an optimal solution has been found. The final tableaux are 





3] 





(6) : 
g 
: 
(7) 2 
3S 





When presented in the form of a Gantt chart, this schedule is 


machine 





$2 





The complete solution tree is illustrated in figure 6. The 
numbers in parentheses beside the tableaux illustrate the correspondence 


between the branches in figure 6 and the tableaux. 


C. ° SUGGESTIONS FOR EFFICIENT COMPUTATION 

The branch and bound algorithm presented clearly yields the optimal 
solution to the n-job m-machine scheduling problem. Its computational 
efficiency, however, has not been examined. An accurate statement re- 
garding its usefulness in solving a real problem, therefore, cannot 
presently be made. Some possible refinements to the algorithm are 
suggested. 

At each lower bound computation the only arithmetic operations 
performed are addition and subtraction. These operations take place 
only when required by the comparisons of job start and finish times. 
The number of additions can possibly be reduced by computing only as 
many Ts *s as necessary. That is, when any T, is greater than a lower 
bound that has been computed, cease calculation on that branch and 
proceed to a branch which appears to offer the best possibilities. In 
other words, proceed to the branch with the smallest lower bound 
Recent ly computed. It may be necessary LOM reLunh co One On ENeSe 
partially completed branches. Storage availability and computation 
expense would determine the tradeoff between storing a partially 
completed tableau and performing all the arithmetic necessary to 
reproduce a tableau. 

Another possible refinement of the method involves the use of the 
graphically computed upper bound of section III. B. This upper bound 


can be used to eliminate entire branches from the branch and bound 
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figure 6 
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solution tree. As the upper bound computed in this manner represents 
the time to complete one feasible schedule, any lower bound discovered 
to be greater than this time may be disregarded. The branch and bound 
procedure clearly will not permit return to a branch with a lower bound 
greater than the optimal schedule. The optimal schedule will not be 
Known, however, until the algorithm has terminated. The use of the 
graphical upper bound can then eliminate the requirement for storing 


any information connected with the branches with excessive lower bounds. 


D. COMPUTATIONAL EXPERIENCE 

A computer program for the branch and bound algorithm starts on 
page 38, but the program employs none of the suggestions for improving 
computational efficiency. Only limited computational experience has 
been gained with this program. Very little can be said regarding the 
efficiency of the algorithm since only four problems have been solved 
with the algorithm. The problem sizes (job x machine), the number of 
iterations to arrive at a solution, the core storage requirements (in 
thousands of bytes), and the computing times (in seconds) are shown 
below. The first problem solved is the example problem discussed 


previously. 


PROBLEM NUMBER SIZE NO. OF ITERATIONS STORAGE COMPUTING TIME 


] 3x4 ZS 82 8.40 
2 3x4 44 86 8.65 
3 4x5 283 120 SSR 
4 4x10 25] 350 Za n0c 


So 





The process times and technological orderings for problems 3 and 4 
were generated randomly. Problem 2 is taken from Brooks and White 


[Ref. 6]. 
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V. CONCLUSION 


Two methods for solving the n-job m-machine scheduling problem have 
been suggested. The graphical method for computing an upper bound on 
the time to process n jobs on m machines is simple to perform but does 
not always yield the optimal solution. The method does, however, yield 
an upper and lower bound on this schedule time. The optimal solution 
has been determined, of course, when these two bounds are the same. 

The branch and bound procedure described, on the other hand, always 
produces the optimal schedule if storage capacity and computation time 
permit. Its computational efficiency, however, must be examined in 


greater detail before it can be considered an efficient algorithm. 
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